Natural growth model of weighted complex networks 
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Abstract 

We propose a natural model of evolving weighted networks in which new links are not necessarily connected 
to new nodes. The model allows a newly added link to connect directly two nodes already present in the 
network. This is plausible in modeling many real-world networks. Such a link is called an inner link, while a 
link connected to a new node is called an outer link. In view of interrelations between inner and outer links, 
we investigate power-laws for the strength, degree and weight distributions of weighted complex networks. This 
model enables us to predict some features of weighted networks such as the worldwide airport network and the 
scientific collaboration network. 



1. Introduction 

Weighted networks are represented by a matrix (toy) 
whose entries are weights attached to links (or edges) 
connecting pairs of nodes (or vertices) i and j. If nodes 
i and j are not connected, then we set u)y = 0. In 
this paper, we will only deal with the symmetric case; 
Wij — Wji. Typical examples of weighted networks are 
the worldwide airport network [2, 4] and the scientific 
collaboration network [4, 8]. In the former, a weight 

is the number of available seats on direct flight con- 
nections between airports i and j, and in the latter a 
weight is the number of papers whose coauthors include 
two authors i and j. 

The strength Si of each node i is defined by [5, 6] 



E ■ 

ier(i) 



where T(i) is the set of neighbors of i and fcj = |r(«)| is 
its degree (the number of its neighbors). The strength 
Si is a generalization of degree and tells us quantity of 
node's significance in the network. In the case of the 
worldwide airport network, the strength provides the 
traffic through i. 

In complex networks the probability distribution 
p(k) that any given node has degree k is fundamen- 
tally important. Many complex networks often exhibit 
a power-law degree distribution 

p(k) oc fc~ 7 with 2 < 7 < 3. 

In weighted complex networks the probability distribu- 
tions p(s) and p(w) of strength s and weight w also play 
a relevant role. In [2, 4] it is shown that weighted com- 
plex networks as above have power-law distributions for 
p(s) and p(w). And these are the subject of the present 
paper, too. 

Most previous approaches [5, 6, 7] to modeling of 
weighted networks are based on the assumption that 
new links always emanate from new nodes. However, 
most real complex networks evolve in aquiring new links 
between existing nodes as well. That is, new links need 



not necessarily emanate from new nodes. In the world- 
wide airport network, for instance, new airline connec- 
tions are frequently created between two existing air- 
ports, independently of the creation of new airports. 
Obviously, for scientific progress, collaboration of active 
scientists may play a more important role than collabo- 
ration with newcomers. 

This paper relaxes this stringent constraint by re- 
suming the weighted network model of [5, 6] combined 
with a recent growth model of [10]. So, as in [10], we 
distinguish two types of a newly added link; an outer 
link and inner link. An outer link is the one emanating 
from a new node, while an inner link connects a pair of 
already existing nodes. 

The strength and degree of each node and the weight 
of each link are analytically solved. After that their 
power-law distributions are derived from them. In par- 
ticular, it is shown that the power-law for strength s; 

p(s) oc s~ 7 with 2 < 7 < 3, 

naturally follows from analysis of the model. The case 
7 — > 2 can be explained from the theory simply by as- 
suming that the ratio of inner links is much higher than 
that of outer ones. 

An analogous argument presented here is applied to 
other growth models in which all new links are always 
assumed to be connected to new nodes. In this way 
more natural growth models for weighted networks can 
be obtained. 



2. Model of weighted networks 

We assume throughout that any existing node i is 
connected to a new outer or inner link at a rate propor- 
tional to 



n(«<) 



(i) 



where the denominator is the sum over all of already 
existing nodes. This is a generalization of preferential 
attachment rule, which is fully discussed in [1, 3, 7, 9]. 
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Based on this attachment rule and [5, 6], the model of 
the evolution for a weighted network is divided into the 
following three parts (i)-(iii). 

(i) It starts with a network having a small number of 
nodes, say No nodes, and several links with the 
same weight wo . Without loss of generality we set 
wo = 1 as in [5, 6]. 

(ii) New nodes and new links constantly enter the net- 
work. Whenever a new node enters the network, 
we advance the time step by one. We assume that 
a new node i enters the network at time t = i 
and has the same amount of strength at first (see 
Eq.(4) below). Moreover, a constant quantity of 
new links, say m links, are added per unit time 
step and each has weight wo = 1 at first. New 
links can be either inner or outer. 

(iii) When a node i is connected by a new link, varia- 
tions of the existing weights are induced across the 
network. That is, the weights Wij of links between 
i and all j g r(i), are modified according to the 
rule 

"Vij 



Wi 



Win + K- 



where k is a constant with k > 0. So in this case 
the strength of node i is modified as Si — > Si+n+1. 

Thus weights are updated and the growth process 
is iterated by introducing new vertices until N nodes 
are obtained. Since the time is measured with respect 
to the number of vertices that enter the network, so mt 
links are added up to t = N — No time steps. 

In order to investigate the evolution of weighted 
complex networks, we make use of the mean-field the- 
ory [1, 5, 6, 7]. Hence all the variables Si 7 ki,Wij, and t 
are treated as continuous ones. 

In the next section, strength Sj(i), degree fcj(i) of 
each node i, and weight Wij(t) for each link are ana- 
lytically solved. After that the power-law probability 
distributions for them are routinely deduced. 

3. Evolution of strength, degree and 
weight 

When a new link is connected to node i, its strength 
Si can increase either if the link connects directly to i, or 
if it connects directly to one j of its neighbors T(i). Each 
added link increases the total strength by an amount 
equal to 2 k + 2, which implies 



^Si(i) *2m(l + K)t. 



(2) 



According to Eq.(l) and (ii), the rate at which a 
node i with strength Si acquires more strength is 



where am is a proportionality constant. The parameter 
a regulates the growth of the network. In particular, 
it determines the numbers of inner and outer links, and 
the growth rates of strength and weight. 

The evolution equation for strength Si is thus given 

by 

~dt 



(2/S +!)*(*) 



(3) 



(2« + 2)i 

Note that Eq.(3) is applied not only when node i is con- 
nected by an outer link, but also when it is connected 
by an inner link. It was only applied to outer links in 
[5, 6], where new links are always connected to a new 
node and a was automatically set as one. 

We must determine the initial condition for Sj in or- 
der for (2) to hold. The solution of Eq.(3) takes the 
form 

/£ X (2k+1) Q /(2k+2) 

where the constant A is the common initial value for 
all Si and fcj, i.e., A = Si(i) = h(i), since the initial 
weight w is unity for each new edge. In order to de- 
termine the constant A, we make use of the continuous 
approximation as follows. 

V- , s . V- ft\ (2k+1W(2k+2) 

E«(*) = A E(i) 

i i 

_ ^(2k+1) Q /(2k+2) y^-(2K+l)q/(2K+2) 

^(2k+1W(2k+2) j\-(2K+l)a/(2 K +2) di 

j.1-(2k+1W(2k+2) 
^(2k+1) Q /(2k+2)_ ' 



(2k+ l)a/(2« + 2) 

_ At 

_ 1 - (2/c + l)a/(2 K + 2) ' 

Hence for (2) we must choose the initial condition for 
each Si and hence that for ki as 

A = Si (i) = fej(i) = (2(k + 1) - (2k + l)a)m. (4) 

This is the number of outer links per unit time step, 
too, because they are connected to a new node. Thus 
the solution of Eq. (3) under the initial condition (4) is 



9i(t) = (2(k + 1) - {2k + l)a)m(r) 



£s (2k+1)o/(2k+2) 



, (5) 



amll(si 



am 



for t > i. 

On the other hand, the remainder among m links, 

B = m - A = (2k + l)(a - l)m, 

gives the number of inner links that connect pairs of 
existing nodes. We will use the ratio 

B (2«+l)(a-l) 
A ~ 2(k+1) - (2K + l)a 
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to indicate the size of inner links. Note that A does not 
become zero, because the case 2(k + 1) — (2k + l)a = 
implies that no properly growing networks can be gen- 
erated. 

Conversely, the parameter a is determined by a given 
ratio B/A as 



a = 



(2 k + 2)B/A+(2k+1) 
(2k+1)(B/A + 1) ' 



Next let us examine the condition for a. Obviously, 
the initial condition (4) must satisfy 

< 2()c+l) - (2k+ l)a < 1, 

and the necessary condition for a is 



1 < a < 



2k + 2 
2k + 1' 



(6) 



The original model of [5, 6] corresponds to the case 
a = 1 or B/A = 0, since inner links are forbidden. In 
contrast, the more a approaches (2k + 2)/ (2k + 1), the 
larger B/A we have, and vice versa. If we have B/A > 1 
as it seems to be the case of many real-world networks, 
then a satisfies the inequalities 

4k + 3 2k + 2 

< a < 

4k + 2 2k + 1' 

and inner links are dominant over outer ones. 

Let us put L(t) — ^\ Si(t). Summing up both sides 
of Eq.(3) then leads us to 

dL _ (2K+l)L(t) 
~dt ~ a (2K + 2)t 

This solution L(t) seems to contradict (2) for large t, 
because L(t) oc £(2k+iW(2k+2) anc [ j ts exponent is less 
than one from (6). In order to get rid of this discrep- 
ancy, using the Dirac delta function S, we should replace 
Eq.(3) by 

ds t _ (2K + l)si(t) 
~dt ~ a (2K + 2)t 

+ (2(« + 1) - (2k + l)a)mS(t - i), 

for t > i. 

Next let us consider the degrees ki(t). They have to 
satisfy 



^h(t) w 2mt. 



(7) 



Using Eq.(l), the rate equation for fcj is similarly given 

by 



dh Si(t) Si(t) 
— = f3m— = fi- 
at 



£. Sj ^(2k + 2)*' 



(8) 



where fi is another constant parameter. Note that fi 
is not equal to two, as expected from (7). In this case 
Eq.(8) should be replaced by 



using the delta function and the initial condition (4) for 
hi . 

We want to determine the parameter fi. Solving 
Eq.(8) from (5) and noting (4), we get 



ki(t) = (2(k + 1) - (2k + l)a)m 

[3 (t\ (2k+1)«/(2k+2) 

(2kTT) 



X (2k + l)a \iJ (2k + l)a J 



(2k +1) 

Hence we have 

h(t) ~ {2(k + 1) - (2k + l)a}m 

i 

(- fi „ (2k+1) Q /(2k+2) 

X I (2k+ l)a ^ \i) 

i 

+ ( 1 -(2kTT)^>}- 
Again using the previous approximation 



EG) 



t\ (2 re +l)a/(2«+2) t 

~ 1 - (2k+ l)a/(2« + 2)' 



and (7), the above relation becomes 
2mt « (2(k + 1) - (2k + l)a)mt 

x( 7 g , (_ L_ + 

I (2k + Da VI- (2k+1)«/(2k + 2) / J 

From this we see that the parameter fi must be taken 
as 

fi = (2k + D"- 2k. 
It follows from (6) that the range of fi is 

1 < fi < 2 

as in [10], and the more fi approaches 2, the larger B/A 
becomes. 

Moreover, the solution fcj(i) can be described by 
means of Si as 



*>(<) - '^iT^ 'M 

(2k + lja 

+ (2(K+l)-(2K+l)a) 



(9) 



2«m 



(2k + l)a' 



Therefore, we see that the relation between strength Sj 
and degree fcj is linear. 

Finally we deal with the evolution of weight using a 
similar argument. The weight Wij increases by the ad- 
dition of a new link connected either to i or to j. So the 
rate equation for weight can be expressed by 



dwi S 
~dT 



am 



/ Si ^Wjj 

VV„s/> Si 



Sj Wij ^ OLKWij (t) 
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with the inital condition w i:( (iy ) = w = 1 , where t^ is 
the time when the link connecting i and j is established. 
Alternatively, it may be rewritten by 



dwij aniVij (t) 



+ u>oS(t - ti ). 



dt (k + l)t 
Solving this equation yields 

+ \ an/ 



W, 



i(t)=( 



U 



(10) 



From (6) it follows that the exponent satisfies the 
inequalities 



k olk 2k 
< < 



K+l ~ K+l 2k + 1 

This means that weight grows faster as B/A becomes 
larger. 

Note that the time tij generally satisfies tij > 
max(i, j) for each j £ r(i), due to the existence of in- 
ner links. Combining (5) and (10) leads to the following 
relation on the time ti 



ier(i) 



n ■ 



-oik/ 



if — 1 



for large t, where A is given by (4). Also from this 
relation it follows that the degree \T(i)\ has to grow ac- 
cording as a or B/A is larger. 

Particularly, in case of k = 0, the degree of node i 
at time t is given by 



|r( i )|«(2- a )m(-)" , 
for large t, as shown in [10]. 



4. Power-law distributions 

From the solutions (5), (9), (10) of s,(i), h(t), Wij(t), 
respectively, it is straightforward to deduce the three 
power-law probability distributions for strength, degree 
and weight. Furthermore, we examine the ranges of 
exponents of power-laws as the fraction of inner links 
changes. 



This implies that, as B/A becomes larger, the exponent 
7 S approaches the constant 2. Since (4k+3)/(2k+1) < 3 
always holds, we obtain the power-law distribution for 
strength s as in Section 1 regardless of k. 

B. Degree 

Similarly from (5) and (9) the degree distribution 
p(k) can be obtained as 



p(k) oc (k — a) 7 , 7 = 1 + 



2k 



(2K + l)a 



where a is the second term of the right hand side of (9); 
a = 2Km{2( K + 1) - (2k + l)a}/((2« + l)a). Thus the 
exponent is the same as strength. So we get the same 
range of the exponent 7; 2 < 7 < (4k + 3)/(2k + 1). 

In particular, as in [5, 6], we need not take the limit 
k — ^ 00 (an unrealistic assumption) in order to explain 
the case when 7^2. The analysis of our theory tells 
us that it is the case when the ratio B/A becomes as 
high as possible. 

If k = 0, then a = and 1 < a < 2. Therefore, we 
have the power-law degree distribution p(k) oc k" 1 with 
2 < 7 < 3; the scale- free distribution [1, 3]. 

C. Weight 

The solution (10) yields the probability distribution 
p(w) of weight; 



p(w) oc w 



7w = 1 + 



K+l 



Also from (6) we get the range for the exponent 7 W as 

2 + ^ <7w <2 + -, 

where k is assumed to be positive. 

An analogous argument can be applied to other 
growth models in [9, 11], for example, where all new 
links are assumed to be connected to new nodes. More- 
over, both models of [9, 11] are based on [5, 6]. In this 
way more natural growth models for weighted networks 
can be obtained from them. 



A. Strength 
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